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Assortative and disassortative mixing investigated using the spectra of graphs 
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We investigate the impact of degree-degree correlations on the spectra of networks. Even though 
density distributions exhibit drastic changes depending on the (dis)assortative mixing and the net¬ 
work architecture, the short range correlations in eigenvalues exhibit universal RMT predictions. 

The long range correlations turn out to be a measure of randomness in (dis)assortative networks. 

The analysis further provides insight in to the origin of high degeneracy at the zero eigenvalue 
displayed by majority of the biological networks. 
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I. INTRODUCTION 

Last two decades have witnessed a rapid advancement 
in the field of complex networks M- The prime idea 
governing this framework is to consider a system made 
of interacting units. To categorize and understand real 
world systems based on interacting units, many mod¬ 
els have been proposed, among which Erdos-Renyi (ER) 
random 0], scale-free (SF) [Bj and small world [§] are 
the most popular ones. Further, degree-degree correla¬ 
tions have also been used as one of the key properties 
of networks characterization 0 3413 and is known to 
confer robustness to biological networks 0 - The ten¬ 
dency of (un)like degree nodes to stick together is termed 
as (dis)assortativity. Various social networks are known 
to be assortative while few of the biological and tech¬ 
nological networks have been reported to be disassorta¬ 
tive [13 18]. Despite its importance for real networks, 
(dis)assortativity does not appear in any of the model 
networks discussed above, and is driven by some other 
mechanism, for example reshuffling algorithm j2()f . While 
spectral behaviour of uncorrelated networks have been 
quite well understood [21|, despite real world systems 
being highly correlated [9], such understanding for the 
correlated networks still needs to be developed. 

Spectral graph theory is an established branch of math¬ 
ematics, and eigenvalues of corresponding adjacency ma¬ 
trices are known as finger prints of the underlying graphs 
[2214251 ]. With recent advancement in the network the¬ 
ory, the spectral graph theory, traditionally used in in¬ 
vestigations of random and regular graphs, got extended 
to studies of graphs motivated by real world systems. 
These spectral studies, apart from presenting bounds for 
extremal eigenvalues highlight their importance by relat¬ 
ing them with the various structural as well as dynam¬ 
ical properties of the networks [26j, [23] ■ The studies of 
networks further reveal a key impact of assortativity on 
the extremal eigenvalues [28}, which has been explored in 
context of disease spreading (29) and diffusion processes 
[30j |. thereby exhibiting the importance of spectral stud¬ 
ies of networks for a more comprehensive understanding 


of complex systems. This paper presents a systematic 
analysis of impact of degree-degree correlations on the 
spectral properties of various networks under the random 
matrix theory (RMT) framework. Since its introduction 
in 1960s, in the context of nuclear spectra, the theory has 
been successfully applied to a wide range of complex sys¬ 
tems ranging from the quantum chaos to galaxy [33, [32| . 
Recently, with a spurt in the activities of network frame¬ 
work, the RMT got its extension in analysis of spectral 
properties of various model networks 1351 34J as well as 
those arising from real world systems [35l [361 ]. 


II. METHODS AND TECHNIQUES 

To quantify the degree-degree correlations of a net¬ 
work, we consider the Pearson (degree-degree) correla¬ 
tion coefficient, given as 01 
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where ji,ki are the degrees of nodes at both the ends of 
the i th connection and M represents the total connections 
in the network. 

The random network of size N and average degree ( k ) 
is constructed using the ER model by connecting each 
pair of nodes with the probability p = ( k)/N Q- These 
networks have assortativity coefficient (r) being close to 
zero or exactly zero. To generate the networks with var¬ 
ious assortativity, we use the reshuffling algorithm [20j . 
In this algorithm, after selecting two pairs of nodes ran¬ 
domly, we sort them according to degree. The highest de¬ 
gree node is then connected to the second highest degree 
node with the reshuffling probability p r , which governs 
the (dis)assortative mixing, i.e. we reconnect a high de¬ 
gree node to a (low) high degree one and low degree node 
to a (high) low degree one. With the probability 1 —p r , 
we rewire them randomly. If new connection resulting 
from this rewiring already exists, it is discarded and the 
previous steps are performed. The process is carried out 
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until a steady value of r is attained. For assortative net¬ 
works, the k degree nodes to form a complete graph with 
the value of r being one, the network should have at least 
{k + 1 + 2 n) nodes, where n can be any integer starting 
from 0. As this condition is not satisfied for all the de¬ 
grees present in the network, the network takes a value 
lesser than one. Similarly, the disassortative network can 
have the value of r less than —1.0. 

We make a further note that at high assortativity val¬ 
ues, all the similar degree nodes being connected among 
themselves form groups [20|. As we decrease the assor¬ 
tativity, the connections within the groups of similar de¬ 
gree nodes decrease and the connections between differ¬ 
ent groups of similar degree nodes increase. For disassor¬ 
tative networks, connections between different groups of 
similar degree nodes exist giving rise to a bipartite-like 
structure [20| . 

The SF networks of size N and average degree ( k ) are 
generated using the Barabasi-Albert algorithm by start¬ 
ing with a completely connected network seed and adding 
new nodes one by one which connect with existing nodes 
using the preferential attachment method [l| . 

The networks are represented in the form of adjacency 
matrix by defining Aij = 1, if i and j nodes are con¬ 
nected otherwise A y = 0. For an undirected and un¬ 
weighted network with N nodes, the adjacency matrix is 
N x N symmetric square matrix entailing all real eigen¬ 
values. We denote the eigenvalues as \,i = 1,2..TV and 
Aj < Ai+i and analyse them under the RMT framework. 
The random matrix studies consider two properties of a 
spectra: (1) global properties such as spectral distribu¬ 
tion of the eigenvalues p{ A), and (2)_local properties such 
as eigenvalue fluctuations around A. In RMT, calcula¬ 
tions of spectral fluctuations are done using the unfolded 
eigenvalues \ = iV(Aj), where N{ A) = J^ . p(X) dX is 
the average integrated eigenvalue density [37|. By using 
these unfolded eigenvalues, nearest neighbour spacings 
are calculated as s* = Aj+i — A ». For symmetric random 
matrices with the mean zero and the variance one, the 
nearest neighbour spacing distribution (NNSD) follows 
GOE statistics given as: 
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which shows a level repulsion at small spacing values 
with an exponential fall for larger spacings indicating 
that nearest neighbour eigenvalues are correlated j37j . 
Whereas the spacing distribution of a matrix whose di¬ 
agonal elements are Gaussian distributed random num¬ 
bers and rest of the elements are zero exhibit Poisson 
statistics (P(s) = exp(—s)) indicating that eigenvalues 
are uncorrelated f37j . 

The intermediate of these two distributions can be 
characterized using the Brody equation [HI: 

P/ 3 (s) = exp (—, (3) 


where A and a are determined by the parameter /? as 

0 +1 


A = (1 + f3)a and a = 


Iwb 


The value of 


Brody parameter lies in the range (0 < j3 < 1). The value 
of (3 being 0, indicates the Poisson distribution, where as 
/3 = 1 corresponds to the GOE distribution. Other values 
of [3 indicates that the distribution lies intermediate to 
these two. 

The NNSD provides a correlation measure of subse¬ 
quent eigenvalues, whereas the A 3 (L) statistic measures 
how the eigenvalues which are L distance apart are cor¬ 
related, and can be estimated using the least-square de¬ 
viation of the spectral staircase function representing av¬ 
erage integrated eigenvalue density N( A) from the best 
fitted straight line for a finite interval of length L of the 
spectrum given by [32J : 


A.s(L; x) = — min f [N(X) — aX — b] 2 dX (4) 

L a,b J x 


where a and b are regression coefficients obtained after 
least square fit. Average over several choices of x gives 
the spectral rigidity, the A 3 (L). For the GOE statistics, 
the A 3 (L) depends on L in the following manner: 


A 3 (P)~4 lni (5) 

For the network spectra considered in this paper, there 
is no analytical form of N, and we perform unfolding 
by numerical polynomial fitting using the smooth part of 
the spectra by discarding eigenvalues towards both the 
ends as well as degenerate eigenvalues, if any 31, 32j. 
This renders the dimension of the unfolded eigenvalues 
less than the dimension of the network. 


III. RESULTS 

The bulk part of the spectra of ER random net¬ 
works with r value being close to zero, follow the well 
known semi-circular law f39l . jioj ] (Fig. [2(f)). The ex¬ 
tremal eigenvalues deviate from the random matrix pre¬ 
dictions and indeed provide various information about 
structural and dynamical properties of corresponding 
systems 0i [ 23 , m [il| . In the following, we present 
results pertaining to the impact of assortativity on the 
spectral properties of networks. It turns out that with 
an increase in the assortativity, the semi-circular distri¬ 
bution, as observed for the uncorrelated ER random net¬ 
works, remains unchanged (Fig. [2(a)-(e)). The largest 
eigenvalue exhibits an increasing trend as already dis¬ 
cussed in [H, 0|- As network is rewired entailing dis- 
assortativity, spectral distribution (/?(A)) acquires a very 
different structure than those of the assortative networks. 
The networks start exhibiting a high degeneracy at zero, 
with overall spectra resembling a double humped struc¬ 
ture (Fig. [2(h)), which becomes more pronounced as the 
disassortativity becomes higher or value of r becomes 
more negative (Fig. [lji)) . This increase in disassorta¬ 
tivity is also accompanied with more number of degener¬ 
ate eigenvalues at zero. There could be various reasons 
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for this high degeneracy, few of them, appropriate in the 
present context are: First, as discussed that disassorta- 
tivity supports bipartite-like structure (Jo) and a com¬ 
plete bipartite network has all zero eigenvalues except 
two. Hence bipartite-like behaviour of the disassorta- 
tive networks presents one of the reasons for the occur¬ 
rence of high degeneracy at zero. Second, tree-like struc¬ 
ture has been demonstrated to yield degeneracy at zero 
eigenvalue 42] and disassortativity encourages tree-like 
structure [20j . which in turn indicates high degeneracy 
at zero. We remark that for large IV, the limiting shape 
of p(A) is known for various cases, which for sufficiently 
dense matrices, tend to follow the Wigner semi-circular 
law typical for the Gaussian matrix ensembles [3§, 0(|, 
whereas an ensemble of sparse random matrices of finite 
size are known to yield states beyond the semi-circular 
law in the tails of the distribution [43l-[45j. For sparse 
random graphs, i.e matrices with 0 and 1 entries having 
smaller p values, while the density distribution p( A) of 
an ensemble exhibit singularities, with the height of the 
peaks being the corresponding multiplicities, the bulk is 
still shown to comply with random matrix predictions 
of Wigner’s semicircular law [ill [47}. Moreover, investi¬ 
gations of various model networks mimicking real world 
properties have revealed that the spectra of these net¬ 
works exhibit degeneracy at zero [33] , as observed for 
the sparse random matrices. On that account, despite 
degeneracy at zero, the bulk of the assortative networks 
following the semi-circular distribution, is not surprising. 

As the spectral density only provides a global be- 
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FIG. 1: (Color Online) Spectral density for Erdos-Renyi ran¬ 
dom networks with different values of assortativity coefficient 
r. All graphs are plotted for the networks with size N = 1000 
and connection probability p = 0.01, averaged over twenty 
different realizations of the networks. 



FIG. 2: (Color Online) The NNSD for Erdos-Renyi random 
networks with different values of assortativity coefficient r. 
All graphs are plotted for the networks with size N = 1000 
and connection probability p = 0.01. Histograms are from the 
data points and solid line is for fitting with Brody distribution 
(Eq.EO). 


haviour of eigenvalues, in order to get insight into lo¬ 
cal fluctuations, we further analyse the short range and 
long range correlations in eigenvalues. The NNSD fol¬ 
lows GOE statistics of RMT (Eq. [2]) for all the values 
of r except for very high values corresponding to the 
highly assortative networks (Fig. [2]). What is interest¬ 
ing that the values of r for which p( A) exhibits a very 
similar behaviour, except a change in the value of the 
largest eigenvalue, the NNSD captures crucial structural 
changes reflected through the value of the Brody param¬ 
eter. For the highest achievable value of the assortativity 
coefficient for the particular network parameter for which 
results are presented, the value of f3 comes out to be close 
to 0.3 (Fig. 0(a)), and as the assortativity decreases we 
witness a smooth transition to the GOE statistics with 
value of the /3 turning one. Depending upon the network 
size, average degree and degree sequences, the highest 
achievable value of r for that network may be different 
(as discussed in the Section II), which might lead to a dif¬ 
ferent value of /?. Fig. H[a)-(d) depict that a very small 
change in the value of r is capable of entailing a pro¬ 
found change in the statistics, in-fact it approaches from 
the Poisson to the GOE. Since a very small randomness is 
known to be enough in introducing the short range corre¬ 
lation in eigenvalues fill ] , for a very small deviation from 
the highest assortativity entails GOE statistics. Since 
the assortativity in network supports a the groups hav¬ 
ing similar degree nodes and as assortativity decreases, 
these distinct groups of nodes observed for very high val¬ 
ues of r gets destroyed leading to a transition from the 
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Poisson to the GOE statistics. As soon as the value of 
r is decreased, and sufficient random connections among 
the groups of similar degree nodes are induced, the value 
of Brody parameter f3 becomes one and no further signa¬ 
ture of structural changes on value of f3 is found with a 
further decrease in the assortativity. 

For disassortative networks which are characterized 
with negative values of r, what is remarkable is that de¬ 
spite these networks displaying distinguishable spectral 
distributions than those of the assortative networks, the 
NNSD yields the value of the Brody parameter (/? = 1) 
bringing them into the universality class of GOE. This is 
not surprising as NNSD is analysed by taking the non¬ 
degenerate part of the spectra, and high degeneracy at 
a particular value, for instance at zero, does not account 
for any effect in the NNSD. As long as the underlying 
network has some random connections, the NNSD dis¬ 
plays the GOE statistics j48j • We remark that all the 
networks considered here form a single connected cluster 
as for disconnected networks, even though each individ¬ 
ual sub-network follows GOE statistics, the spectra taken 
together may lead to a different spacing statistics [HJ. 

In order to get a further insight to the structural 
changes arising due to the changes in r values, we probe 
for the long range correlations in eigenvalues for those 
sets which yield the /3 value one. We find that for all 
these values of r, the long range correlations, measured 
using the A 3 (L) statistic (Eq. 0) follow the universal 
GOE statistics as given by Eq. [5] for a certain value of 
L (denoted as Lq) and deviates from this universality 
afterwards (Fig. [3]). 

Note that a regular network, for instance 1-d lattice 
with a periodic boundary condition, follows Poisson dis¬ 
tribution. As connections are rewired, thereby increasing 
the randomness in the network, value of the Brody pa- 




FIG. 3: (Color Online) The A 3 (L) statistic for Erdos-Renyi 
random networks with different values of assortativity coef¬ 
ficient r. All graphs are plotted for the networks with size 
N = 1000 and connection probability p = 0.01. Solid line is 
the prediction from GOE statistics (Eqs. [4] and 0 and open 
circle are calculated from the network. 


rameter increases with an increase in the rewiring prob¬ 
ability and becomes one at the onset of the small-world 
transition, demonstra ting that nearest neighbour eigen¬ 
values are correlated J48(. For such a small change in 
the network structure there is no visible change in the 
density distribution, but the Brody distribution detects 
even such a small change in the number of random con¬ 
nections, and hence has been proposed to be used as a 
measure of randomness at a fine scale ji^]. After the 
Brody parameter attains a value one, the A 3 (L) statis¬ 
tic has been shown to measure the randomness (in terms 
of Lq in this paper), in the underlying network [49j.As 
rewiring probability increases further, the value of Lq for 
which A 3 (L) statistic follows RMT predictions increases, 
demonstrating that the eigenvalues which are Lq distance 
apart are also correlated. Since Lq provides a measure 
of randomness in a network [49| , for the networks under 
investigation in the present work, it turns out that the 
highest assortative network is least random, as value of 
Lq is least for that particular r value (Fig. 0a)). As as¬ 
sortativity of the network is decreased, the randomness 
of the network increases reflected in the higher value of 
Lq. This increase in the size of Lq continues up to r 
being zero, supporting the fact that network reaches to 
the maximum randomness. The value of Lq then remains 
steady for a further decrease in the value of assortativity 
to the minimum possible value of r, i.e. to the maximum 
disassortativity (Fig. 0c)). As most of the real world 
networks have been reported to posses certain level of 
disassortativity [9), based on the A 3 (L) results we can 
argue that real world systems attempt to have more ran- 
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FIG. 4: (Color Online) (a)-(c) represent the NNSD, (d)-(f) 
present the A 3 (L) statistic and (g)-(i) depict the spectral 
density distribution of ER random networks. All graphs are 
plotted for the networks with size N = 2000, (fc)=10 and for 
average over twenty different realizations of the network. 
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FIG. 5: (Color Online) Spectral density for scale-free net¬ 
works with different values of assortativity coefficient r. All 
graphs are plotted for the networks with size N = 1000, 
(fc}=10, for twenty different realizations. 


domness, thereby leading to be disassortative. What fol¬ 
lows that as value of r increases, by keeping network size 
and average degree same, the value of L 0 for which A 3 (L) 
statistic follows RMT predictions increases, indicating an 
increased amount of randomness in the underlying net¬ 
work. Fig. [4] demonstrates that the behaviour of various 
spectral properties remain unchanged as network size in¬ 
creases. Figs. IU(a)-(c) indicate that the value of Brody 
parameter /? becomes one with a very small decrease in 
the value of r. With a further decrease in the value of r, 
the value of Lq for which the A 3 (L) statistic follow GOE 
statistic increases indicating an increase in the random¬ 
ness as discussed earlier. With a further decrease in the 
value of r in the disassortativity regime, there occurs a 
peak at zero eigenvalue which becomes more pronounced 
as network becomes more dissociative which is also ac¬ 
companied with the deviation from the semi-circular dis¬ 
tribution at very low value of r. 

Further, in order to demonstrate the robustness of the 
universal RMT predictions against changes in the net¬ 
work architecture, we present results for the SF networks 
for various values of r. For r being close to zero, the den¬ 
sity distribution of SF networks exhibit the triangular 
shape [40[ . which, with an increase in the assortativity, 
tends to display flattening of the peak. The range of 
the distribution also shrinks as the assortativity increases 
(Fig. (5l(a)-(e)). On the other hand, as we decrease as¬ 
sortativity, i.e. make the network more disassortative, 
the shape of density distribution starts changing from its 
signature triangular distribution, with the peak at zero 
eigenvalue being more pronounced (Fig. E]T)-(g)). As 
we further increase the disassortativity, the eigenvalues 


PPI networks 

N 

ro 

No(PPI) 

o' 

II 

V. 

£ 

O 

J- 

II 

5- 

i? 

H. pylori 

709 

-0.243 

317 

115 

152 

C. elegans 

2386 

-0.183 

1354 

465 

1124 

S. cerevisiae 

5019 

-0.088 

976 

717 

1149 

H. sapiens 

2138 

-0.084 

864 

423 

643 

D.melanog aster 

7321 

-0.083 

2311 

1389 

1975 

E. coli 

2209 

-0.012 

487 

487 

497 


TABLE I: Comparison of number of zero eigenvalues of PPI 
networks of different species and their corresponding configu¬ 
ration models. ?'o denotes the value of the assortativity coeffi¬ 
cient for the PPI networks. Nq(PPI) denotes the number of 
zero eigenvalues in the spectra of the PPI networks. No(r = 0) 
stands for degeneracy at zero for configuration model with 
r = 0, whereas No(r = ro) denotes the same for the configu¬ 
ration models taking r values equal to the corresponding PPI 
network. 


distribute themselves symmetrically and adopt a double- 
hump shape for highly disassortative networks, clearly 
visible in Fig. W) which is accompanied by a high peak 
at the zero eigenvalue similar to that of the ER random 
networks. It is noteworthy that for highly disassortative 
networks, the spectral density of ER and SF model net¬ 
works behave similarly, deviating from their respective 
signature distributions. Further, the (3 value exhibits a 
transition from the Poisson to the GOE statistics with 
a decrease in the r value. Despite the overall spectral 
density being different from that of the ER networks, 
the NNSD and A 3 (L) statistic display similarity in be¬ 
haviour, which is in line of the argument that the eigen¬ 
values fluctuations are calculated from the smooth ho¬ 
mogeneous part of the spectra by not taking degeneracy 
into account and density is not known to be a real test 
of GOE statistics [50 ]. 

We would like to remark here on the impact and re¬ 
liability of network size considered in the present inves¬ 
tigation. In RMT, different quantities are calculated by 
averaging an ensemble of matrices. However for real sys¬ 
tems, calculations are made as running averages over part 
of the whole spectrum. The random matrix predictions 
can be applied to real world systems if above two are 
equivalent, a property known as ergodicity. More explic¬ 
itly, it means that all members of the ensemble, except 
for a set of measure zero, satisfies the above equivalence 
[5ll 152 ]. Due to the ergodicity, one can construct ma¬ 
trix ensembles in different ways: (a) large dimensional 
random matrices with less number of realizations, or (b) 
smaller dimensional matrices with large number of real¬ 
izations. We consider an ensemble of twenty network re¬ 
alizations with a large dimension, which is already shown 
to be good enough to study various structural properties 
of networks, such as degree distributions, clustering co¬ 
efficients etc. 3. Moreover, individual entities of each 
ensemble follow RMT predictions for NNSD with a good 
accuracy, characterized by % 2 values. As we increase the 
realizations, accuracy increases (Fig. [2] and Appendix). 
Consideration of an ensemble consisting of many more 
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number of network realizations would not lead to sig¬ 
nificant betterment or difference in the following prop¬ 
erties of the network spectra: (1) the Brody parameter 
smoothly turning one with a decrease in the value of r 
at a very fine scale; ( 2 ) a further decrease in the values 
of r leading to an increase in the value of L for which 
spectra follows GOE statistics; and (3) increasing height 
of the peak at zero eigenvalues with an increase in the 
disassortativity, owing to the bipartite-like structure of 
the network. 

Next, in order to investigate if the degree-degree cor¬ 
relations in real world system have different spectral 
behaviour than those of the model networks discussed 
above, we consider the protein-protein interaction (PPI) 
networks of six different species. These networks have 
already been shown to follow universal RMT predictions 
of GOE statistics [53j]. We concentrate here on the oc¬ 
currence of high degeneracy at the zero eigen value. The 
assortativity coefficient and fraction of degenerate eigen¬ 
values are tabulated in Table I. As all the PPI networks 
possess negative value of r as well as have a high de¬ 
generacy at zero, we expect disassortativity to be one of 
the factors governing the degeneracy in the real world 
networks. In order to probe more into the correlation 
between disassortativity and degeneracy at zero, we com¬ 
pare the corresponding configuration model for all PPI 
networks presented above (Table I). It is clearly indi¬ 
cated that as soon as the value of r takes a negative 
value (close to the corresponding PPI network), while 
keeping all other parameters of the system same, there is 
an increase in the degeneracy at zero eigenvalue. 


IV. DISCUSSION AND CONCLUSIONS 

The density distribution of the random networks for r 
value being zero follow the Wigner semi-circular distribu¬ 
tion. Even with change in the assortativity (0 < r < 1), 
the bulk part of the spectra keeps displaying semi-circular 
distribution (Fig.[lja)-(f)), whereas increase in the disas¬ 
sortativity (—1 < r < 0 ) leads to the double hump, which 
is symmetrically distributed around a peak at zero eigen¬ 
value (Fig. [Hi))- The height of the peak increases with 
the increase in the disassortativity of the network. 

The NNSD of the networks with the various 
(dis)assortativity values (1 < r < — 1 ), reveal that there 
is a smooth transition in the /?-value around the very high 
assortativity regime. For very high assortativity values, 
/3 values lie close to zero, and as network becomes less 
assortative /? progresses to one. It might be due to the 
reason that the networks with the highest assortativity 
has groups of similar degree nodes which get perturbed as 
r decreases by making random connections among these 
different groups. For rest of the assortativity values, the 
/3 remains fixed at 1, which corresponds to the universal 
GOE distribution as r value goes to negative end. 

Further, the property of Brody parameter being able 
to detect changes in network structure at a fine scale and 


the increase in Lq of A 3 (A) statistic after attainment 
of /3 value one have several implications, one of them 
that concerns the present work is that the value of (3 
distinguishes two networks based on random connections 
present, while the other is that more assortativity in the 
network corresponds to less randomness. Decreasing the 
assortativity leads to an increase in randomness which 
continues up to the value of r = 0 , for which the network 
is most random ( Lq value being maximum). Then the 
value of L for which A 3 (L) statistic follows GOE predic¬ 
tion starts decreasing and remains steady for a further 
decrease in the value of assortativity up to the minimum 
possible value of r (i.e. up to the maximum disassorta¬ 
tivity case). The SF networks also exhibit the similar 
statistics of eigenvalue fluctuations as for ER random 
networks, where density distribution for r = 0 and for 
lower |r| values show triangular distribution instead of 
semi-circular. Both the networks, however, exhibit high 
degeneracy at zero for the disassortative networks. By 
considering different PPI networks, we further demon¬ 
strate the role of disassortativity governing the appear¬ 
ance of degeneracy at zero eigenvalue. 


In spectral graph theory, most of the works concen¬ 
trate on extremal eigenvalues [22]; whereas the RMT re¬ 
search focuses on distribution of various spectral proper¬ 
ties of random matrices with an extension to the random 
graphs, largely ignoring many graphs properties exist¬ 
ing in real world systems. The analysis carried here is a 
step towards bridging this gap by considering two most 
popular tools of random matrix theory, i.e. density and 
spacing distributions, to understand the impact of one 
of the important properties of graphs i.e. assortativity. 
This property has been increasingly realized as a char¬ 
acteristic of a system puUbfij . Our analysis is another 
demonstration of the importance of spacing analysis in 
understanding impact of degree-degree correlation on the 
network detected through the spectra as for very minute 
changes in r, there are no visible changes in the spectral 
density, but this leads to a very drastic changes in the 
eigenvalues fluctuations demonstrating the impact of r 
values on randomness in a network. 


Furthermore, Aa(L) statistic provides an insight into 
the reason why social networks tend to be assortative 
while biological and technological networks tend to be 
disassortative. As randomness, measured in terms of Lq 
for which the A 3 (A) statistic follows RMT prediction, in¬ 
creases with a decrease in r. A direct implication of this 
result can be witnessed in case of social networks where 
entities are known to be associated in ordered fashion 
(people with similar age or educational profile are often 
more connected) (57|, thus providing a probable reason 
as to why social networks tend to assume an assortative 
topology. On the other hand, it has been reported that 
most of the biological and technological networks possess 
a certain level of disassortativity [7, 3, HH . Also the bi¬ 
ological networks, for instance the PPI networks exhibit 
varying amounts of randomness in their underlying net¬ 
works detected through different values of Lq for which 
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the A 3 (L) statistic follows GOE statistics [ 53 J. This ran¬ 
domness has been attributed to mutations occurring in 
course of evolution [58j]. Relating the disassortative na¬ 
ture of the PPI networks and the randomness they pos¬ 
sess, with the results obtained from our analysis of the 
model networks, suggests that biological networks tend 
to become more disassortative in order to comply with 
their underlying randomness. 

To conclude, we present a systematic analysis of 
the spectral properties of the networks with varying 
(dis)assortativity. We find that assortativity has a pro¬ 
found impact on the spectral properties of the underlying 
networks. At a very high assortativity regime, even with 
a slight decrease in the value of r, the Brody parame¬ 
ter smoothly turns one. A further decrease in the values 
of r leads to an increase in the value of Lq of A 3 (A) 
statistic for which the spectra follows GOE statistics. As 
Brody parameter /3 captures the changes in assortativity 
coefficient at a fine scale [48] and L 0 at large scale (iff , 
which further suggest that when r decreases, random¬ 
ness increases. With a further decrease in r, at around 
r = 0 , the density distribution start exhibiting peak at 
zero eigenvalue which becomes more pronounced as r de¬ 
creases further. Interestingly, most of the studies on net¬ 
work spectra report that the bulk part of the spectra 
of the networks having Gaussian and scale-free degree 
distribution follow semi-circular and triangular distribu¬ 
tions [39!, 4(j respectively, but for highly disassortative 
networks, the spectral density of both the degree distri¬ 
butions can have entirely different behaviour. 

Recently, the realm of assortativity has been realized in 
understanding adaptive synchronization fool , which com¬ 
bined with our results of varying amount of randomness 
for various values of r can be explored further to under¬ 
stand dynamical processes on networks. Further, Table 
I indicates that disassortativity is one of the factors con¬ 
tributing to degeneracy at zero. Prevalence of zero de¬ 
generacy has been implicated in terms of gene duplication 
1591 ]. This, along with the impact of change in topology 
of a network, brought upon by assortativity, leading to 
a profound change in the spectral density, provides a di¬ 
rection to explore the evolutionary origin of real world 
systems ftiol |61], Lastly, since randomness or random 
connections in a network has already been emphasized 
for proper functioning of corresponding systems [621], the 
profound role of assortativity parameter revealed through 
the sophisticated random matrix technique is not only 
important for network community attempting to model 
complex systems, but is interesting for random matrix 
communities at the fundamental level as well. 


V. APPENDIX 


Numerical calculations pertaining to assortative mix¬ 
ing, eigenvalues calculations and A 3 (X) statistic are done 
using FORTRAN code written by the Authors. The 
eigenvalues are calculated by calling LAPACK (Lin- 




FIG. 6: (Color Online) (a), (b) and (c) plot average NNSD 
for an ensemble of twenty realizations for different values of 
r, whereas (a), (b) and (c) plot the ensemble having different 
number of the network realizations. The histogram is drawn 
using the data fro the networks, and the solid line is the fitted 
Brody distribution. For all the graphs N = 1000 and (k) = 
10 . 


ear Algebra PACKage) subroutines into the FORTRAN 
code. The calculation of spacings and polynomial fittings 
are done using MATLAB. 

We present the y 2 values as a measure of goodness of fit 
of the model to data, a lower of y 2 indicating a better fit¬ 
ting. As depicted from Fig. [G] the y 2 values consistently 
decrease with an increase in the number of network real¬ 
izations in the ensemble implicating increase in the accu¬ 
racy reaching to the value of y 2 being less than one lying 
in the acceptable range [63]. For assortative networks, as 
less as three realizations in the individual ensemble are 
good enough to bring y 2 within the acceptable range, 
whereas for r taking negative values, the number of re¬ 
alizations in the ensemble increases little bit more (five 
as depicted in Fig. [Bfc)) in order to bring y 2 within the 
acceptable range. This happens as for disassortative net¬ 
works, there is high degeneracy at zero eigenvalue leading 
to less effective dimension of the unfolded spectra (refer 
discussions in the Methods and Techniques section) and 
hence more number of realizations of the networks are 
required. 
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